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We study (d — l)-dimensional excitations in the d-dimensional color code that are created by 
transversal application of the Rd phase operators on connected subregions of qubits. We find 
that such excitations are superpositions of electric charges and can be characterized by fixed-point 
wavefunctions of (d— l)-dimensional bosonic SPT phases with (Z 2 )®'* symmetry. While these SPT 
excitations are localized on (d— l)-dimensional boundaries, their creation requires operations acting 
on all qubits inside the boundaries, reflecting the non-triviality of emerging SPT wavefunctions. 
Moreover, these SPT-excitations can be physically realized as transparent gapped domain walls 
which exchange excitations in the color code. Namely, in the three-dimensional color code, the 
domain wall, associated with the transversal i ?3 operator, exchanges a magnetic flux and a composite 
of a magnetic flux and loop-like SPT excitation, revealing rich possibilities of boundaries in higher¬ 
dimensional TQFTs. We also find that magnetic fluxes and loop-like SPT excitations exhibit non¬ 
trivial three-loop braiding statistics in three dimensions as a result of the fact that the phase 
operator belongs to the third-level of the Clifford hierarchy. We believe that the connection between 
SPT excitations, fault-tolerant logical gates and gapped domain walls, established in this paper, can 
be generalized to a large class of topological quantum codes and TQFTs. 


I. INTRODUCTION 

Classification of fault-tolerantly implementable logical 
gates in topological quantum error-correcting codes is an 
important stepping stone toward far-reaching goal of uni¬ 
versal quantum computatiorP®. Characterization of log¬ 
ical operators is also essential in understanding braiding 
and fusion rules of an yoni c excitations arising in topolog¬ 
ically ordered system^^. Although some classes of two- 
dimensional topological quantum codes are restricted to 
possess string-like logical operators only, there exist non¬ 
trivial topological quantum codes with two-dimensional 
logical operators. Namely, in the two-dimensional color 
code, transversal membrane-like phase operators lead 
to non-trivial action on the ground space which in duce 
a non-trivial automorphism exchanging anyon labelP^. 
Also, in three or more dimensions, the color code admits 
transversal non-Clifford logical phase gates, which are 
indispe nsab le ingredient of fault-tolerant quantum com- 
putatioii^. 

While characterization of excitations and classification 
of logical operators are intimately related, excitations 
arising from logical phase gates in the color code have 
not been studied. In this paper, we study excitations cre¬ 
ated by transversal phase gates in the color code. Some¬ 
what surprisingly, we hnd that such excitations can be 
characterized by bosonic symmetry-protected topologi¬ 
cal (SPT) phases, which have been actively discussed 
in condensed matter physics communitylSHIIl, Formally, 
the system with SPT order has certain on-site symmetry 
G and its non-degenerate ground state does not break 
any of the symmetries. Studies of SPT phases have pro¬ 
vided better understanding of various quantum phases of 
matter, including topological insulators, the topological 
gauge theory and gauge/gravitational anomalies. How¬ 
ever, studies of SPT phases have yet to find interesting 


applications in mantum information science except for 
few instanceJi®^. 

Our first result concerns an observation that {d — 1 )- 
dimensional excitations, created by transversal Rd phase 
operators in the d-dimensional color code, can be char¬ 
acterized by (d — l)-dimensional bosonic SPT phases 
with (Z2)®'^ symmetry. Namely, by writing the emerging 
wavefunction as a superposition of excited eigenstates, 
we find that its expression is identical to a fixed-point 
wavefunction of a non-trivial SPT phase. The on-site 
symmetries emerge from parity conservation of electric 
charges in the color code. While these SPT excitations 
are localized on (d—l)-dimensional boundaries, their cre¬ 
ations require quantum operations acting on all qubits in¬ 
side the boundaries, reflecting the non-triviality of SPT 
wavefunctions. 

Although these SPT excitations are pseudo-excitations 
(not being eigenstates), they may emerge as gapped do¬ 
main walls in the Heisenberg picture. Our second result 
concerns an observation that application of transversal 
logical gates to a part of the system creates a gapped 
domain wall. In the two-dimensional color code, there 
is a one-to-one correspondence between transparent do¬ 
main walls and two-dimensional logical gates with non¬ 
trivial automorphism of anyons. In the three-dimensional 
color code, gapped domain walls, created by i?3 opera¬ 
tors, transform a loop-like magnetic flux into a composite 
of a magnetic flux and a loop-like SPT excitation, reveal¬ 
ing rich possibilities of boundaries in higher-dimensional 
TQFTs. 

We also study braiding statistics of loop-like excita¬ 
tions in the three-dimensional color code. Our third 
result concerns an observation that loop-like magnetic 
fluxes and loop-like SPT excitations exhibit non-trivial 
three-loop braiding statistics. Namely, if a magnetic flux 
and a loop-like SPT excitation are braided while both of 
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them are pierced through a magnetic flux, the resulting 
statistical phase is non-trivial. The non-trivial three-loop 
braiding statistics results from the fact that the three- 
dimensional color code admits a fault-tolerantly imple- 
mentable logical gate from the third-level of the Clifford 
hierarchy. We also find that excitations, which may con¬ 
dense on the domain wall, exhibit trivial three-loop braid¬ 
ing statistics, implying that domain walls and boundaries 
in three-dimensional TQFTs may be classified by three- 
loop braiding statistics of magnetic fluxes and loop-like 
SPT excitations. 


While the discussion in this paper is limited to a 
very specific model of topological quantum codes, we 
believe that our characterization is more generically ap¬ 
plicable. Namely, we anticipate that in a large class of 
topologically ordered systems, pseudo-excitations result¬ 
ing from fault-tolerantly implementable logical gates can 
be characterized by SPT wavefunctions. We further ex¬ 
pect that these SPT excitations possess non-trivial multi¬ 
excitation braiding statistics and provide useful insight 
into classification of gapped boundaries. We thus view 
results in this paper as a stepping stone toward establish¬ 
ing the connection between characterizations of gapped 
domain walls, fault-tolerant logical gates and braiding 
statistics of SPT excitations. 


This paper is organized as follows. In section [TTj 
we describe string and membrane operators in the two- 
dimensional color code. In section III we show that a 
loop-like excitation in the two-dimensional color code can 
be characterized by an SPT wavefunction with Z2 ® Z2 
symmetry. In section EYl we argue that SPT excitations 
can be physically realized as transparent domain walls 
in the Heisenberg picture. In section we describe 
string, membrane and volume operators in the three- 
dimensional color code. In section lYll we show that a 
membrane-like excitation in the three-dimensional color 
code can be characterized by an SPT wavefunction with 
Z2 ( 8 ) Z2 ( 8 ) ^2 symmetry. In section ivnl we study the 
transparent domain wall in the three-dimensional color 
code. We also study the three-loop braiding statistics of 
magnetic fluxes and loop-like SPT excitations. 


II. MEMBRANE-LIKE OPERATORS IN 
TWO-DIMENSIONAL COLOR CODE 


We begin by considering the two-dimensional color 
code defined on a three-valent and three-colorable lat¬ 
tice A where qubits live on vertices. Colors are denoted 
by A^B,C. An example of such a lattice is a hexago¬ 
nal lattice shown in Fig. where plaquettes are colored 
in A, B, C such that neighboring plaquettes do not have 
the same color. The Hamiltonian is given by 

= (1) 

p p 


fXl (Z) 

where P represents a plaquette, and Sp ^, Sp ^ are ten¬ 
sor products of Pauli-X, Z operators acting on all qubits 
on a plaquette P. Interaction terms Sp , Sp, commute 
with each other for all P,P', and thus the system is a 
stabilizer Hamiltonian. Namely, a ground state \ip} sat¬ 
isfies stabilizer conditions S^j^'^l'ip) = S^p'^\ip) = +\ip) for 
all P. 



FIG. 1: The two-dimensional topological color code. The 
Hamiltonian H is a sum of X-type and .Z-type plaquette 
terms on every plaquette. An open line 7 ^®, consisting of 
thick edges of color AB, defines a string-like operator which 
creates a pair of anyonic excitations on shaded plaquettes of 
color C. 


Anyonic excitations in two-dimensional topologically 
ordered spin systems are characterized by string opera¬ 
tors. To construct them in the color code, we assign color 
labels AB, BC,CA to edges of the lattice A depending 
on color labels of two adjacent plaquettes. Consider a set 
of edges of color AB which form a one-dimensional line 
ryAB (Pig yv"e define 

:= (g) A,, := (g) Z,. (2) 

If is an open line, they commute with all the interac¬ 
tion terms except stabilizers on plaquettes of color C at 
the endpoints of ^ab- Thus, applications of X^^\^ab 
and Z^^\^AB create magnetic fluxes me and electric 
charges ec respectively. The correspondence between 
anyon labels and string operators may be represented as 
follows: 


X^^\^AB me, Z^^\^AB ee- ( 3 ) 


Similarly string operators can be constructed from open 
lines consisting of edges of color BC,CA, 

which lead to 


X^^Ubc ITT'A, ^^*^1730 SA, 


( 4 ) 


There are important subtleties in the above character¬ 
ization of anyonic excitations. First, anyonic excita¬ 
tions with three different color labels are not indepen¬ 
dent from each other since applications of Pauli X, Z 
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operator on a single qubit create composites of anyons 
mAmBrnc^eAe-BS-c respectively. In other words, the fol¬ 
lowing fusion channels exist: 

niA X ruB X me = 1, ba x gb x ec = 1- (5) 

Second, an electric charge ba exhibits the non-trivial 
braiding statistics with tob, but not with This is 

because and always commute with 

each other for any choice of "fBC and j'bc- To fully cap¬ 
ture the braiding statistics in the two-dimensional color 
code, it is convenient to construct an isomorphism be¬ 
tween anyons of the color code and those of the toric 
code. Let ei,mi and 62,7712 be anyons in two decoupled 
(j.e. non-interacting) copies of the toric code. Then the 
following correspondence is an isomorphism which pre¬ 
serves braiding and fusion rules: 

mA ^ 7771, rriB -o- m 2 , ba 62 , Bb ^ ei. ( 6 ) 

In fact, it is known that, on a closed manifold, the two- 
dimensional color code is equivalent to two decoupled 
copies of the toric code under a local unitary transfor- 
matiorPiHll! other words, they belong to the same 
topological phas^^. 

The two-dimensional color code possesses not only 
string-like operators, but also transversal membrane 
(two-dimensional) operators. Let H be a Hadamard op¬ 
erator which exchanges Pauli X and Z operators: 

nxn"^ = z, nzT-O = x. (7) 

The Hamiltonian H is symmetric under transversal con¬ 
jugation by Hadamard operators 

UHU^ =H, H :=<^ Hj. (8) 
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Since % transforms X-type string operators into Z-type 
string operators and vise versa, it exchanges electric 
charges e and magnetic fluxes m\ 

BA ^ mA, GB^mB, mA ^ ba, mB ^ Bb- (9) 

The color code admits another interesting transversal 
menbrane operator. Let us define phase operators, act¬ 
ing on a qubit, by R{9) := diag(l,e*®). Of particular 
importance is the so-called Rm phase operator 

:= diag(l,exp(z7r/2'""^)). (10) 

The i ?2 operator exchanges Pauli X and Y operators: 

R 2 XRI = Y, R 2 YRI = -X. (11) 

Let H be a projector onto the ground state space of the 
color code Hamiltonian. Recall that the lattice A is bi¬ 
partite and qubits can be split into two complementary 
sets T and T‘^. Let us define the following transversal 


(two-dimensional) phase operator 

Th ■■= (S) R2 \j (S) ( 12 ) 

jer jeT'^ 

Then the ground state space is invariant under transver¬ 
sal application of R 2 operators: i? 2 n = ni? 2 ll. This 
two-dimensional membrane operator implements the fol¬ 
lowing exchanges of anyon labels (an automorphism): 

Ba Ba, Bb bb, mA mABA, mB mBBB- 

(13) 

In general, in two-dimensional topologically ordered spin 
systems described by TQFTs, transversal membrane¬ 
like (two-dimensional) operators may induce an automor¬ 
phism of anyon labels which preserves braiding and fusion 
rules (z.e. monoidal center s of c ategories which define 
(2 -I- 1)-dimensional TQFTs)P^^. 

We conclude this section by recalling quantum in¬ 
formation theoretical motivations to study transversal 
membrane-like operators in two-dimensional topologi¬ 
cally ordered spin systems. In quantum information sci¬ 
ence, one hopes to perform quantum information process¬ 
ing tasks in a protected codeword space of some quantum 
error-correcting code. The gapped ground state space of 
topologically ordered systems is an idealistic platform for 
such purposes. But how do we perform quantum compu¬ 
tations inside the protected subspace? Ideally one hopes 
to perform logical operations in a way which does not 
make local errors propagate to other spins. Namely, 
one hopes to perform logical operations by transver¬ 
sal unitary gates acting on each spin as a tensor prod¬ 
uct. Thus, it is important to classify transversally im- 
pleme ntabl e logical gates in quantum error-correcting 
code^^HmH 

Transversal membrane operators, such as 
H and R 2 in the topological color code, are examples of 
fault-tolerantly implementable logical gates as they may 
have non-trivial action on the ground state space (if it 
is degenerate). Our goal is to characterize excitations 
arising from fault-tolerantly implementable logical oper¬ 
ators, and the present paper is dedicated to studies of 
those in the topological color code. 

III. SPT EXCITATIONS IN 
TWO-DIMENSIONAL TOPOLOGICAL COLOR 
CODE 

In this section, we study loop-like excitations created 
by parts of a membrane phase operator R 2 in the two- 
dimensional color code and show that they are character¬ 
ized by a wavefunction of a one-dimensional bosonic SPT 
phase with Z 2 0 Z 2 symmetry. Our finding reveals that 
these loop-like excitations in the color code can be viewed 
as a path integral formulation of an SPT wavefunction, 
leading to a physically insightful proof that such a wave- 
function cannot be prepared by symmetry-protected lo¬ 
cal unitary transformations. We note that the circuit 
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depth of preparing SPT wavefuncti ons w as previously 
studied by using different approache^^^^. 

In this section, for simplicity of discussion, we assume 
that the lattice A is supported on a sphere so that the 
system has a unique ground state IV'gs)- 


A. Loop-like excitation from membrane operator 

To begin, let R{9) be a phase operator R{9) := 
diag(l, e*®). An application of R{9) on a qubit at a vertex 
V creates an excited wavefunction: 


tion basis states. Let "Hno-flux be the fluxless subspace 
of the entire Hilbert space where \tj}) G Rno-Rux satisfies 
'W = +1^) for all P- Let no be the total number 
of plaquettes on the lattice A. We define the excitation 
basis states \pi, ■. ■ ,Pno) G "Hno-flux by 

S^^^\pi,---,Pno) = (l-2pj)|pi,---,P;o) (21) 

for j — 1,..., no- More explicitly, we define them by 

Ipg • • • .Pno) ^ n (1 + ) |0 • • • )• (22) 

i=i 


\m) = m\v • 1 ^ 9 .)- ( 14 ) 

Since R{9) is diagonal in the computational basis, it cre¬ 
ates excitations which are associated with AT-type stabi¬ 
lizers on three neighboring plaquettes Pi,P 2 , P 3 of three 
different colors. We would like to characterize this wave- 
function in the excitation basis: 

1 ^ 0 )) ^ ^Pl,P2,P3 \Pl,P2,P3) (15) 

P1.P2.P3 

where Pi,p 2 ,P 3 = 0,1 and |pi,P 2 ,P 3 ) represents an eigen¬ 
state of the Hamiltonian H with 

-5'pfV'liP2,P~3) = (1 - 2pi)|p'i,P'2,P'3) 
S^]^^\pi,P2,P3) = (1 - 2p2)b'l,p'2,p'3) (16) 

S^]^^\P1,P2,P3) = (1 - ‘^P3)\Pl,P2,P3)- 

In other worlds, |pi,P 2 ,P 3 ) records the presence or ab¬ 
sence of electric charges at Pi, P 2 , P 3 while there is no 
other excitations in the system. One hnds 

|V'(6')) h^cos(6»/2)|0,0,6) -kfsin(6»/2)|i,i,i). (17) 

So the phase operator R{9) corresponds to 

R{9)^ R{9)=exp(^i^Xp,Xp,Xp^'^ (18) 

where Xp^^Xp^^Xp,^ act on |0) and |1) as Pauli X op¬ 
erators. Here we used “ 1 —)•” to denote the map from real 
physical systems to the excitation basis. 

Now consider a subset of qubits V and a restriction 
of the membrane phase operator i ?2 onto V, denoted by 

R 2 \v- 

(8) (g) (i?2li)-'. (19) 

j&vnT j&vnT'= 

Consider an excited wavefunction 

li’v) ~R2\v\i’gs)- ( 20 ) 

We hope to represent \tpv) in the excitation basis. First 
let us be more precise with a definition of the excita¬ 


Importantly, not all the basis states are physically al¬ 
lowed since there are certain Z 2 constraints on values of 
Pj. Observe that 

n s?’ = n s?’=n=<8) 

PeA PeB PeC Vj 

where A, B, C represent sets of plaquettes of color A, B, C 
respectively. Let Nji, N 13 , Nc be the total number of 
excitations on plaquettes of color A, B, C respectively. 
Then, 


Nji = Nis = Nc (mod 2). 


(24) 


In other words, an excitation basis state \pi,...,Pno) 
is physically allowed if and only if 'Zij-Pj^APj ~ 
T.j-.p^eBPj = T.j:P^ecPj modulo 2. (Indeed, if 
\pi,---,Pno) does not satisfy this condition, then the 
righthand side of Eq. (22) becomes zero). Because the 
system is not degenerate, basis states |pi,... ,Pno) with 
Eq. (241 span the fluxless subspace Hno-flux completely. 


Thus, the excited wavefunction can be characterize as 
follows 


\tpv) ^ E ^pi,...,PnQ Itl ; ■ • ■ 5 Ptiq) 


(25) 


Pl,...,Pn 


in the excitation basis where is a complex num¬ 

ber with proper normalization and pi,..., pn^ satisfy 

Eq. pll. 


Let us then study a loop-like excitation created by 
a part of the membrane-like i ?2 operator. Consider a 
set of plaquettes of color C which form a contractable 
connected region of qubits V with a single boundary 
(see Eig. [^. Consider an excited wavefunction = 
R 2 \v ■ Since i? 2 |f/’gs) = If/'gs)? the phase operator 

R 2 \v creates excitations only around the boundary of V. 
A key observation is that i? 2 |v creates excitations only 
on plaquettes of color A, B because V can be viewed as 
a set of plaquettes of color C. Namely, let dV be a set of 
plaquettes of color A, B on the boundary (Fig. [^. Ob¬ 
serve that, on the boundary dV, plaquettes of color A 
and B appear in an alternating way. We denote bound¬ 
ary plaquettes by Ai, Bi,..., A„, where 2 n is the to¬ 
tal number of boundary plaquettes. Then the excitation 
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wavefunction can be written as 

IV'W^l4'y)®|0---6) (26) 

where ) = J2p>^p\pm,PBi,---,Pa„,PbJ- Here, the 
first part represents excitations on the boundary 

and the second part |0,..., 0) represents the rest. 



FIG. 2: A loop-like excitation created by R 2 phase operators. 
Filled dots represent qubits in V and plaquettes crossed by a 
closed loop around V form dV. The region V is constructed 
from a set of plaquettes of color C. One applies R 2 operators 
on filled circles and on filled double circles. 

The next task is to find an expression of |44). From 
Eq. the excited wavefunction is given by 

IV'y)^ n exp(±z^M,lfc) • |6)®"“ (27) 

where ± in the product corresponds to i ?2 and (i? 2 )~^ 
in the phase operator i ?2 respectively. Here {i,j,k) rep¬ 
resents a vertex shared by three neighboring plaquettes 
i,j,k. After careful calculations, one can find that the 
boundary wavefunction, in the excitation basis, is given 

by 

l4y) = ^2|oy|6)®2”. (28) 

where 

n 

U2\dv JJexp ■ exp . 

i=i 

(29) 

The boundary wavefunction \4'^gy) can be viewed as 
a one-dimensional system of 2 n qubits supported on a 


closed loop. It is worth finding the Hamiltonian which 
has the boundary wavefunction \ (I>^qy) as a unique ground 
state. Recall that Zp10)®"“ = |6)®"° for all plaquettes P. 
Interaction terms in the Hamiltonian are then obtained 
by considering U 2 \avZpU 2 \gy One finds 

Hov = - '^Xb^.^Za^Xey 

3 3 

(30) 

One can verify that is the unique gapped ground 

state of the Hamiltonian Hgy. This Hamiltonian for 
the boundary wavefunction is identical to that of the so- 
called cluster state up to transversal application of the 
Hadamard operators to each and every qubit along the 
boundary: 

2n 

FIcluster = - 51 ZyiX,Z, + i = (3^) 

i=i 

The ground state (the cluster state) is specified by 

Z,yX,Z,+i\^} = IV') for all j. (32) 

In quantum information science community, the cluster 
state is known as an important resource state for realizing 
the measurement-based quantum computation scheme^. 


B. SPT excitation with Z2 (8)^2 symmetry 

The cluster state is perhaps the simplest example of 
one-dimensional bosonic SPT phases with Z 2 G) ^2 sym¬ 
metry. Let us define a pair of Z 2 on-site symmetry oper¬ 
ators as follows: 

n n 

5a := (g) X2j Sb ■■= (g) X2,_1 . (33) 

i=i 

By on-site, we mean that symmetry operators are 
transversal. The cluster state I'i/’duster) is the unique 
ground state of the Hamiltonian (Eq. ( [^ ) and is sym¬ 
metric under Sa and Sb- 

5a It^cluster) — 5 b |'?/:cluster) — IV^cluster) ■ (^4) 

This can be verified by noticing Y[j=even Zj-iXjZjpi = 
Sa and and 

^j—l^j^j-\-l\'4^cluster} — ['^cluster)- Lct |'0trivial) — 

be a trivial product state with Z 2 0 Z 2 symmetry where 
|-|-) := 4 ( 1 *^^ symmetric wavefunctions 

IV^ciuster) and IV'triviai) are connected by a local unitary 
transformation and belong to the same quantum phase 
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in the absence of symmetries. Indeed, one has 

2n 

l^cluster) ~ exp ^ — It^trivial)- 

i=i 

On the other hand, in the presence of Z 2 0 Z 2 symmetry, 
they belong to different SPT phases. Namely, there is no 
symmetry-protected local unitary transformation U such 

that [ff, tS. 4 ] = [t/, tSi?] = 0 and jV^cluster) — fflV^trivial)- 

For instance, one sees that the unitary transformation 
in Eq. ( p^ does not commute with Sa or Sb- In this 
sense, a cluster state is a non-trivial SPT wavefunction 
with Z 2 0 Z 2 symmetry. 

It is interesting to observe that loop-like excitations 
in the two-dimensional color code are characterized by 
wavefunctions of a one-dimensional SPT phase. To un¬ 
derstand what this observation means, let us further es¬ 
tablish the connection between SPT phases and loop-like 
excitations. Indeed, on-site symmetry of the boundary 
wavefunction naturally emerges from parity con¬ 

servation on the number of anyonic excitations in the 
topological color code. Since R 2 phase operators are ap¬ 
plied on qubits supported on plaquettes of color C, there 
is no excitation associated with plaquettes of color C. 
This implies 


A/a = A/b = 0 (mod 2) (36) 


which leads to the following on-site symmetries in the 
excitation basis 


5a =n =n ■ 


(37) 


After transversal Hadamard transformation, these sym¬ 
metry operators are identical to those in Eq. (33). As 


such, excitations arising in the two-dimensional color 
code are natural platforms for constructing wavefunc¬ 
tions with Z 2 0 Z 2 symmetry. 



FIG. 3: An SPT excitation localized on the boundary dV. 
Creation of this loop-like excitation on dV requires operators 
acting on all qubits in V. If a magnetic flux m crosses an SPT 
excitation, it gets transformed into a composite of a magnetic 
flux and an electric charge e. 

The remaining question then is why the boundary 


wavefunction corresponds to a non-trivial SPT phase. 
The key observation is that, while these loop-like excita¬ 
tions are localized along the boundary dV, they cannot 
be created by a local unitary transformation acting on 
physical qubits in the neighborhood of dV. To see this, 
let us create a pair of magnetic fluxes which are located 
outside of V and move one of them inside E, crossing 
the SPT excitation (Fig. [^. Since i ?2 phase operators 
exchange Pauli X and Y operators, the magnetic flux m 
will be transformed into a composite of an electric charge 
e and a magnetic flux m upon crossing the SPT excita¬ 
tion. Let € be a string operator corresponding to the 
propagation of a magnetic flux into V in the absence of 
an SPT excitation. (It is a tensor product of Pauli X 
operators). Suppose that there exists a local unitary U 
which creates an SPT excitation by acting only on qubits 
in the neighborhood of dV. Then UiW differs from i 
only at the intersection with the boundary. This implies 
that a magnetic flux remains to be a magnetic flux inside 
V, leading to a contradiction. Thus, to create a loop¬ 
like SPT excitation, one needs to apply a local unitary 
transformation on all qubits inside V (or all qubits in the 
complement oiV). 

This argument enables us to show that two symmet¬ 
ric wavefunctions IV^ciuster) and I'i/'triviai) belong to differ¬ 
ent SPT phases. This is because symmetry-protected 
local unitary operators in the excitation basis have 
some corresponding local unitary operators in the two- 
dimensional color code. Suppose that there exists a 
symmetry-protected local unitary transformation U such 
that It/’ciuster) = Cf IV'triviai), which Can be gcnerically Writ¬ 
ten as 


tj = r 



ex.p{—iHsp(t))dt 


(38) 


where [Hsp{t),SA] = [Hsp{t),SB] = 0 and Hsp{t) is geo¬ 
metrically local and consists only of terms with bounded 
norms. Here T represents the time-ordering. Note Hsp{t) 
can be written as a sum of symmetry-protected local 
terms such as Xa^Xa^, , Za^, Zb^ and their 

products. These symmetry-protected local unitary op¬ 
erators correspond to some local operators in the topo¬ 
logical color code. Namely, Xa^Xa^ and Xb^Xb^ corre¬ 
spond to string operators of color BC and CA, consisting 
of Pauli Z operators, which end at plaquettes Ai, Aj and 
Bi,Bj respectively. Moreover, Za^ and Zb^ correspond 
to X-type plaquette operators of color A, B. As such, 
the existence of a symmetry-protected local unitary U 
implies the existence of a local unitary transformation U 
which acts only on physical qubits near the boundary, 
yet creates a loop-like SPT excitation. This leads to a 
contradiction. Thus, we can conclude that IV'ciuster) ^^nd 
IV^triviai) belong to different SPT phases. 

We have argued that an excited wavefunction \(j)^gy) 
corresponds to a non-trivial SPT phase due to the Z 20 Z 2 
parity conservation and the non-trivial automorphism of 
anyons induced by the R 2 phase operator. In general. 
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loop-like excitations may involve plaquettes of three dif¬ 
ferent colors when V is chosen to be an arbitrary con¬ 
nected region of qubits. The excited wavefunction then 
can be viewed as a one-dimensional system with spins of 
three different colors which possesses Z 2 0 Z 2 symmetry 
with respect to SaSb SaSc where Sa,Sb,Sc are ten¬ 
sor products of Pauli Z operators acting on spins of color 
A,B,C in the excitation basis. As such, our characteri¬ 
zation of excitations by SPT phases is valid in these cases 
too. 

For readers who are familiar with the literature of SPT 
phases, it may be clear that the characterization of a 
wavefunction in the excitation basis in the fluxless sub¬ 
space "Hno-flux is essentially equivalent to “ungauging” the 
on-site symmetries, opposite to the procedure of gaug- 
ing the on-site symmetries of SPT wavefunctionsI^Eoi^ 
In this picture of ungauging, the two-dimensional color 
code with a loop-like excitation serves as a path inte¬ 
gral formulation of a one-dimensional SPT wavefunc¬ 
tion. Namely, creation of a loop-like SPT excitation via 
transversal R 2 operators on the bulk V can be interpreted 
as a symmetry-protected quantum circuit preparing a 
non-trivial SPT wavefunction with Z 2 0 Z 2 symmetry. 
An interesting application of this picture is that, if one 
considers the two-dimensional color code embedded on 
a hyperbolic surface, one obtains a MERA (multi-scale 
entanglement renormalization ansatz) circuit for an SPT 
wavefunction with Z 2 Z) 

Finally we remark that similar SPT excitations emerge 
in two decoupled copies of the toric code where the 
transversal control-.^ operator preserves the ground state 
space. In this setting, a loop-like SPT excitation involves 
electric charges from two copies of the toric code, each 
possessing one copy of the Z 2 symmetry. This conclusion 
also follows from the unitary equivalence of the color code 
and two decoupled copies of the toric code on a closed 
manifolcP^. 


IV. GAPPED DOMAIN WALLS AND 
FAULT-TOLERANT LOGICAL GATES 

While one-dimensional SPT excitations, created by i ?2 
phase operators in the two-dimensional color code, are in¬ 
teresting from a theoretical viewpoint, they do not exist 
as stable objects since they are superpositions of eigen¬ 
state excitations which would decohere immediately. In 
this section, we argue that SPT excitations are physi¬ 
cally realized in a certain way. Namely, we point out that 
SPT excitations can be viewed as transparent gapped do¬ 
main walls in the color code in the Heisenberg picture. 
Our finding also reveals an intriguing relation between 
classifications of domain walls and fault-tolerantly im- 
plementable logical gates. 

We note that boundaries in the two-dimensional quan¬ 
tum double model are discussed in Ref.l^. We also note 
that domain walls in SPT phases are discussed in Ref.l^. 


L domain wall R 



FIG. 4: Construction of a transparent gapped domain wall in 
the topological color code. R 2 operators are applied only to 
the qubits on the right side of the system (a shared region). 
Anyons get transformed when crossing the domain wall as 
depicted in the figure. 

A. Transparent domain wall 

The key idea is to transform the Hamiltonian by 
transversal operators instead of transforming the ground 
state. To begin, consider the two-dimensional color code 
supported on the honeycomb lattice (Fig. |^. Recall that 
the transversal Hadamard operator, 77 = pre¬ 

serves the Hamiltonian. Here we split the entire system 
into two parts, the left part L and the right part R. Con¬ 
sider the restriction of the transversal Hadamard opera¬ 
tor onto the right part of the system R: 

HU = (g)77|. (39) 

v^R 

and the transformed Hamiltonian H = ('H|_R)i7('H|_R)^. 
Note that the resulting Hamiltonian H remains the same 
as before except on the boundary between L and R: 

H = Hl + Hr + Hlr (40) 

where Hr and Hr are the same as the corresponding 
terms in H. Since H-Ir is a unitary transformation, the 
transformed Hamiltonian H remains gapped. Here, Hrr 
can be viewed as a transparent gapped domain wall con¬ 
necting Hr and Hr, which swaps the electric charge e 
and magnetic flux m as follows: 

(e^|m^), (toaIba), (eslms), (msles). (41) 

By a transparent wall, we mean that no single anyon from 
either side of the wall may condense on the wall. Here, 
(a|5) represents that an anyon a gets transformed into b 
by crossing the wall from the left to the right. Trans¬ 
parency of the wall imposes that (a|l) or (l|o) for a ^ 1 
is not allowed where 1 denotes the vacuum. Similarly, 







the transversal R 2 phase operator on R creates a gapped 
domain wall which changes anyon labels as follows: 

(toaIgawa), (caIca), (mslesmB), (esles). (42) 

This transparent domain wall corresponds to the loop¬ 
like SPT excitation created by R 2 operators in the 
Heisenberg picture as shown in Fig. The ground state 
space of the topological color code is also symmetric un¬ 
der the so-called T transformation: 

T : X^Y, Y ^ Z, Z ^ X. (43) 

Transversal application of T operators leads to the fol¬ 
lowing transparent domain wall: 

(caItoa), (toaIgawa), (gb|tob), (mslGsTOs)- (44) 


Clearly, the aforementioned construction of transpar¬ 
ent gapped domain walls works for arbitrary transversal 
membrane-like operators in (2 -|- l)-dimensional TQFTs 
as long as the operators induce non-trivial automorphism 
among anyon labels. It turns out that the construc¬ 
tion works not only for transversal ope rator s but also for 
any locality-preserving transformation^^^. Formally, a 
locality-preserving unitary transformation U is defined 
to satisfy the following condition. For an arbitrary uni¬ 
tary operator £ supported on some region H, there ex¬ 
ists £' supported on H U dV such that \UIU^ — £'| « 0 
where dV is the boundary of V of finite width. Examples 
of locality-preserving transformations include local uni¬ 
tary transformations with removing and adding ancilla 
qubits. Also translating all qubits on the lattice by finite 
sites is a locality-preserving transformation. By applying 
such a transformation with non-trivial automorphism of 
anyon labels on the half of the system, one can create a 
transparent gapped domain wall which changes anyon la¬ 
bels according to the automorphism associated with the 
membrane-like operator. 


To give a concrete yet non-trivial example of locality¬ 
preserving transformation, consider the two-dimensional 
toric code supported on a square lattice on a torus: 
Htoric = — ~ Ep ^p- Observe that transversal ap¬ 

plication of Hadamard operators, followed by shifting all 
the lattice sites in a diagonal direction, leaves the Hamil¬ 
tonian invariant. This transformation is clearly locality¬ 
preserving and swaps electric charge e and magnetic flux 
m. If one applies this transformation partially on one 
side of the system, one is able to create a domain wall 
which exchanges e and m upon crossing the wall. (Equiv¬ 
alently, one may consider the Wen’s formulation of the 
toric code where unit translations lead to exchange of e 
and rrP^. We note that this domain wall in the toric 
code was previously constructed by BombiiJ^. 


B. Membrane operator and domain wall 


We have seen that membrane operators associated with 
non-trivial automorphisms among anyon labels lead to 
transparent gapped domain walls in (2 -|- l)-dimensional 
TQFTs. An interesting question is whether the presence 
of transparent gapped domain walls implies membrane¬ 
like locality-preserving transformations. In this subsec¬ 
tion, we present complete classifications of transparent 
gapped domain walls for the two-dimensional toric code 
and the two-dimensional color code. Namely, we find 
that every transparent domain wall has a corresponding 
membrane operator associated with non-trivial automor¬ 
phism of anyon labels. 

In (2-1-1)-dimensional TQFTs, labels of anyons, which 
may condense on a gapped boundary, can be character¬ 
ized by maximal sets of anyons with trivial self and mu¬ 
tual braiding statisticJi^. To characterize a gapped do¬ 
main wall which connects two topologically ordered sys¬ 
tems, we fold the system along the domain wall and view 
the domain wall as a gapped boundary where two sys¬ 
tems (the left and the right) are attached togethei!^. We 
are then able to classify automorphisms of anyon labels 
on the domain wall as condensation of anyons from the 
left and the right in the folded geometry. 

To begin, let us characterize transparent domain walls 
of the toric code. Anyons are denoted by e, m. One finds 
the following two transparent domain walls: 

Wo : (G|e), (m|m) 

Wi : (e|m), (to|g). ^ ’ 

Here Wq represents a trivial domain wall where anyon 
labels are not altered and Wi corresponds to the afore¬ 
mentioned domain wall in the toric code which exchanges 
e and m. One can view the domain wall Wi as a gapped 
boundary which absorbs eirrir and meer where the sub¬ 
scripts i and r denote the left and the right. One then 
sees that and meCr have trivial self and mutual 

braiding statistics. We note that this classification of 
transparent domain walls in the toric code is well known 
in the literature, se^^ for instance. 

Next, let us characterize transparent domain walls for 
the two-dimensional color code, which is unitarily equiv¬ 
alent to two copies of the toric code. The problem can 
be reduced to finding 4x4 matrices with binary en¬ 
tries which satisfy certain conditions reflecting triviality 
of braiding statistics in the folded geometrji^. We found 
72 different types of transparent domain walls in the color 
code by an exhaustive search. Observe that transparent 
domain walls form a group since a product of two walls 
is also a wall. As such, domain walls can be constructed 
from a certain complete set of generators of domain walls. 
Below, we list five types of domain walls which form a 
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complete set. 

Wi : {eA\mA), {mA\eAmA), (eslms), (TOslesms) 

W2 ■■ (e^lmA), (toaIba), (bbItob), (wsles). 

(46) 

These domain walls preserve the color labels of anyons 
while exchanging electric charges and magnetic fluxes. 
Here Wi corresponds to T operators which permute Pauli 
X, Y, Z operators and W 2 corresponds to Hadamard op¬ 
erators. Here Wi and W 2 form a subgroup which is iso¬ 
morphic to the symmetric group S 3 . Thus the group of 
domain walls is non-abelian. The following domain walls 
exchange color labels of anyons: 

W3 : (eAles), (toaIuib), (esIeA), (tobIwa) 

W 4 : (baIbb), (wAlmB), (bbIbaBb), (wBlmAmB)- 

(47) 

One can construct the corresponding membrane opera¬ 
tors by decoupling the color code into two copies of the 
toric code and exchanging them. Finally, we find the 
following domain wall 

W5 : (baItob), (mA|mA), (bbIbb), (tob|ba) (48) 

which corresponds to a domain wall in a single copy of 
the toric code. We find that domain walls Wi,..., W 5 
generate all the 72 different transparent gapped domain 
walls in the color code. The main conclusion is that, for 
every transparent domain wall in the color code, there 
exists a corresponding membrane-like locality-preserving 
operator which preserves the ground state space and in¬ 
duces non-trivial automorphism of anyon labels. 

Classification of gapped domain walls is an impor¬ 
tant problem in condensed matter physics community 
as many of realistic physical systems have boundaries 
and domain walls. Classification of fault-tolerantly im- 
plementable logical gates is of relevance to the quantum 
information science community as they are indispens¬ 
able building blocks for fault-tolerant quantum computa¬ 
tion. Whether the correspondence between domain walls 
and membrane operators generically holds for (2 -|- 1)- 
dimensional TQFTs is an interesting future problem to 
study. 


V. VOLUME OPERATOR IN 
THREE-DIMENSIONAL COLOR CODE 


Consider the three-dimensional color code defined on a 
four-valent and four-colorable lattice A where qubits live 
on vertices. Colors are denoted by A, B, C, D, and are 
associated with volumes. The Hamiltonian is given by 

( 49 ) 

P G 


where P represents a plaquette and G represents a vol- 
(jc) (z) 

ume (Fig. ^|. Here Sq ,Sp ' commute with each other 
due to the four-valence and four-colorability of the lat¬ 
tice. For simplicity of discussion, we assume that the lat¬ 
tice A is supported on (homomorphic to) a three-sphere 
so that the ground state is unique. A systematic proce¬ 
dure of constructing such a lattice is knowiP^. 



FIG. 5: Stabilizer operators in the three-dimensional color 
code. X-type stabilizers are associated with volumes while 
Z-types stabilizers are associated with plaquettes. 


Three-dimensional topologically ordered spin systems 
may possess both point-like and loop-like excitations 
which are characterized by string and two-dimensional 
membrane operators respectively. We begin by con¬ 
structing string operators. Given a four-valent and 
four-colorable lattice A, one can assign color labels 
ABC, ABD, ACD, BCD to its edges since, for a given 
edge, there always exist three volumes of different col¬ 
ors sharing the edge. Consider a set of edges of color 
ABC which form a one-dimensional line such that 

a string connects volumes of color D as shown in 

Fig.ga). We define 


:= (^ Zj. (50) 

j^^ABC 

If is an open line, Z^^^\.yABc commutes with all 

(X) 

the interaction terms except stabilizers Sq on volumes 
of color D at the endpoints of 7abc- Thus, one can 
characterize electric charges as follows 

Z^\^ABc eo, ^ABD ^ ec 

CB, ZBCD\ .yBCD CA- 

To construct membrane operators, we assign color la¬ 
bels AB, AC, AD, BC, BD, CD to plaquettes of the lat¬ 
tice A. Consider a set of plaquettes of color AB which 
form a two-dimensional sheet (membrane) (3^^ as shown 
in Fig.[6];b). We define 

^1/3^- := (g) Xp (52) 

j^pAB 

If is an open membrane with boundaries, the oper¬ 
ator creates excitations associated with stabilizers S^p ^ 
on plaquettes of color CD on boundaries of Pab- Thus, 


























10 


one can characterize loop-like magnetic fluxes as follows 

"mcD, mBD 

X^^\pAo mBC, X^^\pBcruAD ( 53 ) 

X^^\pBn mAC, X^^\pCB VriAB- 

These excitations with different color labels are not inde¬ 
pendent from each other since the following fusion chan¬ 
nels exist 

ca X cb X ec X e^i = 1, biab x mco = 1 

(54) 

ruAC X ruBD = 1, ruAD x uibc = 1- 

It is convenient to construct an isomorphism between 
anyons of the color code and those of the three- 
dimensional toric code. Let ei, toi, 62 ,TO 2 , 63 ,m 3 be 
anyons in three decoupled copies of the toric code. Then 
one has 

ca ^ Cl, cb o 62 , ec 63 

(55) 

TTiAB ^ m3, ruAc ^ m2, mBC ^ mi. 

It is known that, on a closed manifold, the three- 
dimensional color code is equivalent to three decoupled 
copies of the toric code under a local unitary transforma- 
tioipMl. 




FIG. 6: (a) An open line 7 '^^'^, consisting of edges of color 
ABC, which defines a string-like operator Z^^^\..^abc ■ The 
string connects volumes of color D, and creates excitations 
on two volumes of color D sitting at the endpoints, (b) An 
open sheet (membrane) , consisting of plaquettes of color 
AB, which defines a membrane-like operator X^^I^ab. The 
membrane connects plaquettes of color CD, forming a loop¬ 
like flux, and creates excitations on plaquettes of color CD 
which are on the boundary of the sheet . 

Note that the three-dimensional color code has volume 
and membrane phase operators. Recall that the lattice 
A is bipartite and qubits can be split into two sets T and 
For a phase operator R 3 = diag(l, let us de¬ 

fine the following transversal volume (three-dimensional) 
phase operator 

Th-=^R^\o (56) 

jeT jer<‘ 

It is known that R 3 preserves the ground state space of 


the three-dimensional color code. Namely, if there is only 
a single ground state, one has 

Rsllpgs) OC llpgs). (57) 

This is a rather non-trivial statement, and readers are 
referred t<P for a proof. 

We then define membrane-like i ?2 phase operators. Let 
us consider a set of plaquettes of color AB which form a 
two-dimensional sheet . Recall that the sheet is 
bipartite and qubits can be split into two sets fl T 
and [3^^ We define 

:= (g) R2\j (g) {R2\j)-^- (58) 

iG/3'^snr= 

If is an open membrane with boundaries, the op¬ 
erator commutes with all the stabilizers except X-type 

(X) 

stabilizers Sq ’ of color C and D on boundaries. So, 
R2^’'\p AB crCcitiGS cL loop“lik 6 Gxcitcitioii involving both, co* 
and Cd- This membrane-like phase operator R 2 ^^\pAB is 
closely related to the other membrane operator 
by the following relation: 

K(R3,'X^\pab) oc R2^^\0ab (59) 

where K{U,V) := UVWV"^ is the so-called group com¬ 
mutator. One can see that the excited wavefunction, re¬ 
sulting from R 2 ^^\i 3 ab , corresponds to a one-dimensional 
SPT wavefunction with Z 2 0 Z 2 symmetry in the exci¬ 
tation basis where Z 2 O Z 2 symmetries are associated 
with parity conservations of ec and cd- We shall de¬ 
note such a loop-like SPT excitation by scd- Loop-like 
SPT excitations, created by R 2 membrane operators, are 
characterized by 

R 2 ^^\pAB Scd, R 2 ^^\fiAc sbd 

R2^^\i3^d SBC, R2^^\i3BC SAD (5^) 

r) BD \ T) CD \ 

H 2 \^bd ^ saci ^2 l^co ^ Sab 

with the following fusion channels 

Sab X Scd = 1, sac x sbd = 1, sad x sbc = 1- (61) 

For instance, if sab and scd are located next to each 
other, then one can eliminate them by local unitary trans¬ 
formations. 


VI. SPT EXCITATIONS WITH Zj (g) Za (g) Za 
SYMMETRY 

In this section, we study membrane-like excitations 
created by applying the phase operator R 3 in some con¬ 
nected region of qubits. Namely, we show that membrane 
excitations are characterized by a wavefunction of a two- 
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dimensional bosonic SPT phase with Z20Z20Z2 symme¬ 
try. We also find a new simple fixed-point Hamiltonian 
for SPT phases with Z 2 0 Z 2 0 Z 2 symmetry. We conjec¬ 
ture that the corresponding SPT phase is the so-called 
type-III phase which is dual to the non-abelian quantum 
double model based on a dihedral group D^. 


A. Membrane excitations 

Consider a set of volumes of color D which forms a 
connected contractible region of qubits V with a sin¬ 
gle boundary. We consider a restriction of the three- 
dimensional phase operator R-^ on V, denoted by R^\v'- 

i?3|y:= (g) i?3|, (g) (i?3|,)-'. (62) 

j&vnT jevnT'= 

Consider an excited wavefunction \ipv) ■= Rslv ' IV'gs)- 
Since the region V is constructed from a set of volumes 
of color D, Itpv) involves excitations only on volumes of 
color A, B, C. Since the phase operator R^ is diagonal in 
the computational basis, it creates excitations which are 
associated with A-type stabilizers on volumes. We hope 
to represent IV'v) in the excitation basis in a way similar 
to the two-dimensional case. Not all the basis states are 
physically allowed since 

n =n =n =n 

p&A PeB Pec PeT> 

where A,B,C,'D represent sets of volumes of color 
A,B,C,D respectively. Let iVe, Nc, Nx> be the to¬ 
tal number of excitations on volumes of color A, B, C, D. 
Since A-p = 0, one has 


Na = Nb=Nc = 0 (mod 2) (64) 

which leads to Z 2 0 Z 2 0 Z 2 symmetry of j'i/'y) in the 
excitation basis. Clearly, the wavefunction I'lpv) has ex¬ 
citations only on the boundary dV of V since R^ltpgs) oc 
\ijjgs)- So, one can characterize the excitation wavefunc¬ 
tion by the boundary wavefunction \4>^Qy) which corre¬ 
sponds to the presence and absence of excitations asso¬ 
ciated with volumes on the boundary dV of V. 

Let us find the boundary wavefunction \4i^Qy). A key 
observation is that the boundary dV can be viewed as 
a three-valent and three-colorable lattice (see Fig. [^. 
Namely, since F is a set of volumes of D, the boundary 
dV consists of plaquettes of color AD,BD,CD. Con¬ 
sider volumes of color A, B, C which are located on the 
boundary dV. Such volumes on the boundary, which 
form the basis states of \(j)^gy), can be associated with 
plaquettes of color AD,BD,CD on the boundary dV. 
As such, the boundary wavefunction \4>^gy^) can be in¬ 
terpreted as a two-dimensional wavefunction which in¬ 
volve qubits supported on plaquettes of a three-valent 
and three-colorable lattice dV. Here, it is convenient 



FIG. 7: The boundary of a volume of color D, viewed 
as a three-colorable lattice consisting of plaquettes of 
color AD, BD, CD. The boundary excitation wavefunc¬ 
tion \cj>QY) is supported on qubits associated with three- 
colorable vertices on the dual lattice dV*. The fig¬ 
ures were created using Robert Webbs Stella software 
(http: / / www.software3d.com/Stella.php). 


to introduce a dual lattice picture (Fig. [^. Let us con¬ 
struct a dual boundary lattice dV* by viewing the centers 
of plaquettes in dV as vertices. Then, in a dual lattice 
dV*, qubits live on vertices and vertices are colored in 
A,B,C. In summary, the boundary wavefunction {(ji^gy) 
can be viewed as a two-dimensional system of qubits sup¬ 
ported on a triangular lattice without boundaries whose 
vertices are colored in A,B,C (Fig. [^. 

One can find that the boundary wavefunction \(j)^gy) 
is given by 

\4v) = n |0)®” (65) 


where n represents the number of vertices and (i, j, k) 
represents a triangle in the dual lattice dV*. Here “±” in 
the product corresponds to i ?3 and (i? 3 )“^ in the phase 
operator i ?3 respectively. One can construct a Hamil¬ 
tonian which has this wavefunction \4'^gy) as a unique 
ground state: 

H = - exp{iTTtgtq^)^ (66) 

3 {]qq') 

where the product runs over triangles {jqq') containing 
the vertex j , and 


^ ■“ 2 


(67) 


Here Xg is a Pauli-A operator acting on a qubit at a 
vertex q. The boundary wavefunction is a unique gapped 
ground state of this Hamiltonian (assuming dV is a closed 
manifold): 


Qj I (l>av ) = +14 ^^) (68) 
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and has Z 2 ® Z 2 0 Z 2 symmetry emerging from parity 
conservation in Eq. (64). 




FIG. 8: The Hamiltonian for an SPT phase with Z 2 0^2 0^2 
symmetry, supported on a triangular lattice with vertices of 
three different colors. 


B. SPT phase with Z2 0 Z2 < 8 ) Z2 symmetry 


In this subsection, we further establish the connec¬ 
tion between the boundary wavefunction and SPT phases 
with Z 2 0 Z 2 0 Z 2 symmetry. T o be consistent with 
the notation used in the literatur^SE^, we shall apply 
transversal Hadamard operators to the wavefunction and 
the Hamiltonian in Eq. (65) and Eq. (66). After exchang¬ 
ing X and 0, the wavefunction is given by (Fig.[^ 


IV’Spt) := exp I ± — ZiZjZk j |-|-)®” 


(69) 


and the Hamiltonian is given by 


Hspt - X! II *^^99' 

3 <j99'> 


(70) 


where the so-called control-Z operator CZg^/ is given by 

Z 1 

CZqq, := exp{iTTSqSq') Sq := -. (71) 

(Here CZqq> applies a Pauli Z operator on q' provided the 
qubit at g is in a state |l)q). The wavefunction satisfies 
the stabilizer conditions: 


^ilV'SPx) =+|V'spt) Vj. (72) 

Let us define symmetry operators as follows: 

5^ := (g) Sb := (g) Xf \ 5c := (g ^ 

j j j 

(73) 

which are associated with parity conservation. One can 
verify that the wavefunction is symmetric: 

5a|V’spt) = 5b|V'spt) = 5c|V’spt) = IV'spt)- (74) 


To see this, notice that the following relations hold 


i[qj=sa, (75) 

JGA j&B j&C 

which can be shown by recalling the fact that (CZ^g/)^ = 
I. Thus, the boundary wavefunction has Z 2 0 Z 2 0 Z 2 
symmetry. 

Readers who are familiar with the literature on bosonic 
SPT phases may notice that the aforementioned Hamil¬ 
tonian is identical to the fixed-point Hamiltonian of two- 
di mensio nal SPT phases with Z 2 symmetry proposed in 
RefEMl]^ instead of Z 2 0 Z 2 0 Z 2 symmetry. The crucial 
observation is that this Hamiltonian can be also viewed 
as a fixed-point Hamiltonian of two-dimensional SPT 
phases with Z 2 0 Z 2 0 Z 2 symmetry if an appropriate 
set of symmetry operators is considered. In Ref.“ the 
Hamiltonian in Eq. (70), supported on a triangular lat¬ 


tice, was treated as a Z 2 symmetric model with respect 
to SaSbSc- Yet, if one views the same Hamiltonian with 
Z 2 0 Z 2 0 Z 2 symmetry with respect to Sa 0 5^ 0 5c, 
the system can be viewed as an example of SPT phases 
with Z 2 0 Z 2 0 Z 2 symmetry. This change of symmetry 
operators has a highly non-trivial consequence since the 
resulting Z 20 Z 20 Z 2 symmetric wavefunction is different 
from three decoupled copies of Z 2 symmetric wavefunc- 
tions. 

We then argue that the wavefunction IV'spt) belongs to 
a non-trivial SPT phase with Z 20 Z 20 Z 2 symmetry. We 
present two different arguments to reach this conclusion. 
First, one can consider a thought experiment of creating 
a magnetic flux and making it cross this excitation as in 
two-dimensional case. As discussed further in the next 
section, a magnetic flux becomes a composite of mag¬ 
netic flux and a loop-like SPT excitation upon crossing 
the membrane excitation. This implies that membrane 
excitations associated with R 3 operators cannot be lo¬ 
cally created by applying local unitary transformation 
along the boundary. The second argument replies on the 
fact that IV'spt) belongs to a non-trivial SPT with Z 2 
symmetry, which has been already established in Ref.l^ 
by gauging the on-site symmetries. This implies that 
IV'Spt) must be non-trivial with respect to Z 2 0 Z 2 0 Z 2 
too. It remains to determine which class of SPT phases 
IV^SPt) corresponds to. 

It is known that there are 128 different SPT p hases 
with Z 2 X Z 2 X Z 2 symmetry including the trivial on j^^ l ^^ l 
These 128 different phases can be fully characterized by 
three different types of SPT phases called type-I, type-H 
and type-HI. Their corresponding cocycles are given by 


w/(A, B, C) = exp -b c« - [&(*) -b c^*)]) 

0 Jii{A,B,C) =exp ('yaW(5(*) -bcW - [&« -b c^]) 


B, C) = exp 


(76) 
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where l<i<j<£<3. Note A = (a^^\ B = 

( 5 ( 1 )^ ^( 2 )^ ^( 3 ))^ C = and = 

0,1. Here brackets represent modulo 2 calculus. By gaug¬ 
ing on-site symmetries, these three types of SPT phases 
are mapped to the twisted quantum double models which 
are spin systems with intrinsic topological order (with¬ 
out symmetries). The topological model associated with 
type-III is known to be dual to the (untwisted) quantum 
double model based on the dihedral group D 4 with non- 
abelian topological order. We conjecture that the above 
Hamiltonian with Sa(^Sb®Sc symmetry corresponds to 
the type-III class since the wavefunction mixes three dif¬ 
ferent modes A, B, C and possesses three-party entangle¬ 
ment among H, B, C. Another argument supporting this 
conjecture relies on the consideration of boundaries. As 
discussed in a forthcoming papeil^, in the d-dimensional 
quantum double model based on a finite group G, gapped 
boundaries can be constructed by considering d-cocycles 
with respect to H C G. For three decoupled copies of the 
three-dimensional toric code {i.e. the three-dimensional 
quantum double model with Z 2 0 Z 2 0 Z 2 symmetry), 
one can construct a transparent domain wall by using 
3-cocyle of Z 2 ® Z 2 ® Z 2 C (Z 2 0 Z 2 0 Z 2 )®^ which cor¬ 
responds to the type-III class. One can verify that such 
a transparent domain wall transforms excitations as de¬ 
scribed in Eq. ([^ and Eq. ( |78| . It is also worth noticing 
that operators X and R 2 form a projective representa¬ 
tion of D 4 . We note that one may also identify the cor¬ 
responding SPT phase by considering symmetry twist or 
flux insertion as proposed in Ref.HlHSl. 

To conclude this section, we briefly mention excitations 
in the higher-dimensional color code. There exists a d- 
dimensional generalization of the topological color code 
with point-like electric charge and codimension -1 mag¬ 
netic flux. The code admits the transversal Rk phase 
operators (for 1 < k < d) which are supported on k- 
dimensional regions and preserve the ground state space. 
One can then study codimension-1 excitations created by 
d-dimensional Rd-type operators, which lead to a (d— 1)- 
dimensional boundary wavefunction with (Z 2 )®‘^ sym¬ 
metry. The emerging model is supported on a (d — 1)- 
dimensional lattice where qubits are associated with ver¬ 
tices which are d-colorable. The Hamiltonian for this 
wavefunction involves a Pauli X operator and generalized 
(d — l)-qubit control-Z operator^, and the colorability 
of the lattice is essential to ensure that the system has 
(Z 2 )®‘^ symmetry. 


parent domain walls in (3-|-l)-dimensional TQFTs cannot 
be classified by transpositions of eigenstate excitations 
such as electric charges and magnetic fluxes. In this sub¬ 
section, we demonstrate that in the three-dimensional 
color code a transparent domain wall may transform a 
loop-like magnetic flux into a composite of a magnetic 
flux and a loop-like SPT excitation. 

As in the two-dimensional case, we split the entire sys¬ 
tem into two parts, left L and right R. We then apply 
i ?3 operators only on qubits in R. This creates a gapped 
domain wall on the boundary between L and R (Fig. |^. 
Since the i ?3 phase operator commutes with Pauli-Z op¬ 
erators, electric charges cross the domain wall without 
being affected: 


(e^leA), (esles), (ec|ec)- (77) 

Yet, magnetic fluxes get transformed in an interesting 
way. Recall that R^XR\ oc Ai? 2 - Then, one finds that a 
magnetic flux becomes a composite of a magnetic flux and 
a one-dimensional SPT excitation with 'L 2®'^2 symmetry. 
Namely, one has 

[mAB\mABSAB), {'mBc\mBCSBc), imcA\mcASCA) 
(sab|sab), {sbc\sbc), [sca\sca)- 

(78) 

Thus, in order to characterize gapped domain walls in 
the three-dimensional color code, SPT excitations need 
to be considered. 



FIG. 9: (a) A transparent domain wall created by R 3 phase 
operators applied to the right part R of the system, (b) Two- 
loop braiding process, (c) Three-loop braiding process. 


VII. TWO-DIMENSIONAL BOUNDARY AND 
THREE-LOOP BRAIDING 

A. The R3 boundary 

In (2-|-l)-dimensional TQFTs, transparent gapped do¬ 
main walls are classified by automorphisms of anyonic ex¬ 
citations. The most important distinction between two- 
dimensional and three-dimensional systems is that trans- 


B. Braiding statistics between two excitations 

In (2-|-1)-dimensional TQFTs, anyons which may con¬ 
dense on the boundary can be classified by cons idering 
maximal sets of mutually bosonic anyon^l ^Fspe ^ 
argument can be generalized to higher-dimensional sys¬ 
tems too, leading to a conclusion that excitations which 
may condense on boundaries must have mutually trivial 
braiding statistics. Thus, in order to classify boundaries 
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and domain walls in (3+l)-dimensional TQFTs, braiding 
statistics among particles and loops need to be studied. 
In this subsection and the next, we study braiding pro¬ 
cesses among electric charges, magnetic fluxes and SPT 
excitations in the three-dimensional color code. This sub¬ 
section is devoted to braiding statistics of two excitations 
(Fig. I^b)) while the next subsection is devoted to three- 
loop braiding statistics (Fig.j^c)). 

In three-dimensional systems, there are three possi¬ 
ble braiding processes which involve two excitations: (i) 
particle-particle, (ii) particle-loop and (iii) loop-loop. By 
braiding, we mean a process of creating excitations, wind 
one around the other, and annihilating them separately. 
Here we are interested in the resulting 1/(1) statistical 
phase. First, consider braiding between electric charges. 
Since electric charges 6 ^, 65 , ec are mutually bosonic, 
one has 


^i0(eA,eB) _ 


(79) 


where 0 ( 6 ^, &b) represents the corresponding 1 /( 1 ) statis¬ 
tical phases. The above equation also holds under per¬ 
mutations of color labels. 

Second, consider braiding of a particle and a loop-like 
excitation. Braiding statistics is given by 

^i 9 (eA,mBc) — _ ^ie(eA,mcA) _ 


So, electric charge ca and a magnetic flux tubc have 
non-trivial braiding statistics with additional —1 phase 
factor. This is because propagations of are character¬ 
ized by a string-like Pauli-Z operator ^bco while 

propagations of mBC are characterized by a membrane¬ 
like Pauli X operator X^^pAn. Namely, a string 
and a membrane overlap with each other odd times. 
The braiding of SPT excitations and an electric charges 
is trivial since SPT excitations are superpositions of elec¬ 
tric charges: 


gi9{eA,SAB) _ ^i0{eA,SBc) _ ^ 


(81) 


Third, consider braiding of two loop-like excitations. 
In general, braiding statistics of two loop-like excitations 
a and /3 can be probed by a process depicted in Fig.[Io|(a). 
Let Ua be a unitary operator corresponding to creating 
a pair of loops a and a', and sending a' in a direction 
perpendicular to the loop a as depicted in Fig. [^b). 
Let Up be a unitary operator corresponding to creating 
a pair of loops /3 and /3', and enlarge /3' as depicted in 
Fig. 10'c). Note that Ul^ and Up are unitary operators 


whicnreverse the processes of Ua and Up. The following 
sequence of unitary transformations implements a braid¬ 
ing of two loops a and /? (Fig. [Io|(a)): 




(82) 


Recall that K(U^, Up) = U^UpUaUp. As such, two-loop 
braiding statistics of a and p can be characterized by 


the vacuum expectation value of the group commutator 
K{Ul,Up). Note that continuous deformations of the 

braiding trajectories do not change the value of 
and thus braiding statistics are topologically invariant. 

Let us now consider a braiding of a magnetic flux 
and an SPT excitation sbc- Let <5 be a loop, where the 
trajectories of niAB and sbc intersect (see Fig. [Ioj(a)). 
Note that Ua for a = tuab is a membrane operator with 
X operators acting on plaquettes of color CD. Similarly, 
Up for a = SBC is a membrane operator with R 2 op¬ 
erators acting on plaquettes of color AD. As such, the 
commutator K(U^, Up) is a loop 6 of Pauli Z operators 
acting on edges of color ACD. This operator acts as an 
identity operator in the ground state space. Thus, we 
conclude that 

^i0{mAB,SBc) — 2 (83j 

Similar argument allows us to conclude that = 1 

for all loops a, /3. 



FIG. 10: (a) Two-loop braiding process as a group commuta¬ 
tor K{Ul, Up), (b) Pair of loop excitations a, a'. (c) Pair of 
loop excitations P,l3'. 


C. Gapped domain wall and three-loop braiding 

Recently, Wang and Levin have proposed that full 
characterization of three-dimensional SPT phases re¬ 
quires analyzing the three-loop braiding statistic^^. In 
this subsection, motivated by this pioneering work, we 
study three-loop braiding processes among loop-like mag¬ 
netic fluxes and SPT excitations and show that they ex¬ 
hibit non-trivial three-loop braiding statistics. 

The three-loop braiding is depicted in Fig. j^c) where 
one first creates a loop 7 , and then creates/braids a and 
/3, and then makes the system return to the vacuum. 
The key distinction from two-loop braiding is that, a 
and /3 are braided in the presence of another loop 7 . 
We shall denote the corresponding statistical phase by 
giS{a,p,j) ^ The braiding statistic of three loops a,j3,"/ can 
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be probed by the process depicted in Fig. 11 Let U. 


be a unitary operator corresponding to creating a pair 
of loops 7 and 7 ' from the vacuum and enlarging 7 '. 
Then, the following sequence of unitary transformations 
implements the three-loop braiding: 




where 


(85) 


As such, the three-loop braiding statistics corresponds 
to the vacuum expectation value of the sequential group 
commutator. 



FIG. 11: Three-loop braiding process as a sequential group 
commutator K{K{U^, Up), U^). 


Let us now consider the braiding of uiab, sbc, mcA- 
Recall that K{U\,Up) is a string of Z operators, char¬ 
acterized by edges of color ACD. Recall that is a 
membrane of X operators acting on plaquettes of color 
BD. These two operators overlap odd times, and thus 
the group commutator between K{U^,Up)''' and leads 
to a non-trivial statistical phase: 

pie(mAB,SBC,'mcA) — 


In general, one finds 

^ie{sK ,mi(, ,m^„) _ ^ie{mK,SK' ,mK") — _ 

where {K, K', K”} = {AB, BC,CA}, and permutations 
of K, K', K”. For any other combinations of loo-like ex¬ 
citations a, , 0 , 7 , one has = 1 . 

Finally, we study the braiding statistics of particles 
and loop-like excitations in a transparent domain wall 
of the three-dimensional color code in a folded geome¬ 
try. Let us begin with the braiding with two excitations. 
The only non-trivial braiding processes are the ones in¬ 
volving e and m, such as the braiding of (e^je^) and 
{'mBc\'<TiBCSBc)- Braiding of electric charges ca and 
magnetic fluxes ruBC on each side of the wall contributes 
to —1 phases in the braiding statistics, which cancel with 
each other. For the three-loop braiding, the non-trivial 
processes are the ones involving a = {mBc\'^BCSBc), 
b = {mcA\'mcAScA), c = (mAB\mABSAB)- Due to the 
symmetry under permutations of color labels, it suffices 


to study 6 {a,a,a), 6 (a,a,b), 6 (a,b,a), 6 {a,b,c). It is 
easy to verify = I. As 

for 6 {a,b,c), there are two “— 1 ” contributions to the 
phase, and thus one obtains = 1. In conclu¬ 

sion, on a transparent R 3 wall, particles and loops that 
may condense on the wall possess mutually trivial braid¬ 
ing statistics, in two-excitation and three-loop braiding 
processes. This result verifies the hypothesis that gapped 
domain walls and boundaries can be classified by trivi¬ 
ality of braiding statistics, including three-loop braiding 
processes. 


To conclude the paper, we would like to further discuss 
a relation between braiding processes and group com¬ 
mutators. We have seen that two-loop and three-loop 
braiding processes are characterized by a group com¬ 
mutator K{U^,Up) and a sequential group commutator 

K{K{Ul, ,U)^). In the d-dimensional color code, one 
can define codimension -2 excitations with non-trivial d- 
excitation braiding statistics. This process can be char¬ 
acterized by the following sequential group commutator: 


K{...K{Vc„K{Vo,,,Vc,J)...) (88) 

where V^_ corresponds to a unitary operator creating 
excitations of type aj for j = 1,..., d. Readers from 
quantum information science community may find this 
expression particularly interesting since the sequential 
group commutator is essential in defining the Clifford 
hierarchy, the central object in the classification of fault- 
tolerant logical gates in topological stabilizer code^^^^. 
Indeed, the non-trivial three-loop braiding statistics in 
the three-dimensional color code is a direct consequence 
of the fact that transversal i ?3 operator belongs to the 
third-level of the Clifford hierarchy. These observations 
may hint a deeper connection between classification of 
fault-tolerant logical gates, braiding statistics in higher¬ 
dimensional TQFTs and quantum error-correcting codes. 
It will be also interesting to study logical gates and 
boundaries i n top ological quantum codes which are be¬ 
yond TQFTpEa. 
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